ADA  08  77  3  3 


> 


PLASMA  INSTABILITY  IN  ELECTRON  AND  POSITRON 
COLLIDING  BEAMS  IN  STORAGE  RING 


Han  S .  Uhm 

Naval  Surface  Weapons  Center 
White  Oak,  Silver  Spring,  Maryland  20910 

and 


Chuan  S.  Liu 

Department  of  Physics  and  Astronomy 
University  of  Maryland 
College  Park,  Maryland  20742 


4 


aJ 

o 

tu' 


Physics  Publication  Number  80-180 


j  This  document  hns  boon  apr»t>ved  ; 
!  tor  public  rel'»asv  end  sale;  i  s  j 

|  distribution  is  unlimited.  ^ _ j 


UNIVERSITY  OF  MARYLAND 
DEPARTMENT  OF  PHYSICS  AND  ASTRONOMY 

COLLEGE  PARK,  MARYLAND 

80  6  20  160 


li)  :>6 


Plasma  Preprint  PL  #80-040 


(() 


,  _PLASMA  ^STABILITY  IN  JLECTRON  AND^OSITRON 
^  COLLIDING  J&EAMS  IN  STORAGE  RING  f 


W 


ruhm  rA.  o  >  > .  &  //  i  u 

Naval  Surface- Vfcepows  Center  '  _ 

White  Oak,  Silver  Spring,  Maryland  20910 


)>  / 

LHan  S./u 

face  Vfca 


and 


Chuan  S.  Liu 

Department  of  Physics  and  Astronomy 
University  of  Maryland 

College  Park,  Maryland  20742.  - 

A1  QOO/V  ~  77- C  - 


!*! 

Physics  Publication  Number  80-180 
Technical  jfcepaStt^JJumber  80-103 

April  1980 

If )  -  ;  -  '.Pf'i '  />rr-  ■  s- '  /? 


A 


i -Hi  ? 


PLASMA  INSTABILITY  IN  ELECTRON  AND  POSITRON 


COLLIDING  BEAMS  IN  STORAGE  RING 


Han  S .  Uhm 

Naval  Surface  Weapons  Center 
White  Oak,  Silver  Spring,  Maryland  20910 

and 


Chuan  S.  Liu 

Department  of  Physics  and  Astronomy 
University  of  Maryland 
College  Park,  Maryland  20742 


ABSTRACT 


The  f ilamentation  instability  of  the  electron  and  positron 
colliding  beams  in  a  storage  ring  are  investigated  within  the 
framework  of  the  rigid  beam  model  and  the  Vlasov-Maxwell  equa¬ 
tions,  and  closed  algebraic  dispersion  relations  for  the  complex 
eigenf requency  w  are  obtained.  It  is  shown  that  the  typical 
growth  rate  of  instability  is  a  substantial  fraction  of  the 
electron  plasma  frequency  thereby  severely  limiting  the 

electron  density  in  a  storage  ring.  Moreover,  the  influence  of 
collective  self-field  effects  on  the  electron  and  positron  col¬ 
liding  beams  in  the  storage  ring  is  investigated.  The  analysis 
is  carried  out,  distinguishing  the  cases,  where  (a)  the  particle 
motions  are  in  a  very  coherent  orbit,  and  (b)  the  randomness 
dominates  the  operational  condition  of  a  storage  ring  (e.g.,  the 
incoherent  collision  location  by  small  fluctuation,  etc.)  In  — 


2 


I 

I 


INTRODUCTION 


There  is  a  growing  interest  in  the  equilibrium  and  stability  properties 
of  the  electron-positron  colliding  beams  in  a  storage-ring  facility.  A 
recent  experiment^  with  colliding  electron-positron  beams  at  DESY  has  shown 
the  broadening  of  the  beam  cross  section,  thereby  leading  to  reduction  of 
luminosity.  To  address  this  serious  problem,  we  examine  the  f ilamentation 
instability^  of  electron-positron  beams  and  the  influence  of  the  collective 
self-field^  on  the  electron-positron  colliding  beams  in  the  storage  ring. 

For  the  analytic  simplicity,  we  assume  that  beams  have  cylindrical  shape 
and  are  azimuthally  symmetric  in  the  equilibrium  state.  Equilibrium  and 
stability  properties  of  planar  geometric  beans  are  to  be  presentd  in  a 
subsequent  publication. 

O 

In  Sec.  II,  we  treat  the  f ilamentation  instability  of  colliding 
electron-positron  beams  with  finite-geometry  effects  included.  Stability 
analysis  of  dipole  oscillation  is  carried  out  in  Sec.  II. A,  within  the 
framework  of  a  rigid  beam  model,  which  provides  a  simple  instructive 
description.  In  Sec.  II. B,  the  analysis  for  the  high  harmonic  perturba¬ 
tions  with  iK 2  (where  &  is  azimuthal  harmonic  number)  is  carried  out  within 
the  framework  of  the  Vlasov-Maxwell  equations.  An  important  conclusion  of 
the  present  analysis  is  that  the  typical  growth  rate  of  the  f ilamentation 
Instability  is  of  the  order  of  the  electron  plasma  frequency  wpe»  thereby 
severely  limmiting  the  electron  density  in  a  storage  ring.  However,  the 
analysis  of  broadening  of  beam  cross  section  by  repeating  interaction  bet¬ 
ween  electron  and  positron  beams  is  not  completed  yet. 

The  influence  of  the  collective  self-fields^  on  the  electron  and  posi¬ 
tron  colliding  beams  in  the  storage  ring  is  investigated  in  Sec.  III.  The 


theoretical  analysis  is  carried  out,  distinguishing  the  two  cases,  where 
(a)  the  particle  motions  are  in  a  very  coherent  orbit  and  (b)  the  random¬ 
ness  dominates  the  operational  condition  of  storage  ring  (e.g.,  incoherent 
collision  location  by  fluctuation,  etc.).  In  either  case,  it  has  been 
found  that  the  self-fields  effects  play  a  dominant  role  in  the  stability 
behavior  of  transverse  orbit  and  the  expansion  of  beam  cross  section. 
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II.  FILAMENTAT ION  INSTABILITY 


Electron  and  positron  colliding  beams  in  storage  ring  are  likely  sub- 
Ject  to  various  macro-  and  micro-instabilities.  *  Perhaps  one  of  the  most 
important  instabilities  of  the  electron  and  positron  colliding  beam  in  a 
storage  ring  is  the  f ilamentation  instability.  The  unstable  modes 
propogates  nearly  perpendicular  to  the  beam  with  mixed  electrostatic  and 
electromagnetic  components,  the  latter  destabilizing  and  the  former 
stabilizing.  The  perturbed  magnetic  field  is  mostly  in  the  plane  perpen¬ 
dicular  to  the  beam  and  the  Lorentz  force  causes  the  beam  to  filamentate, 
similar  to  the  Weibel  instability.  Unlike  the  Weibel  modes,  which  are 
purely  electromagnetic  for  counter-streaming  electron  beams,  the  linear 
perturbations  of  colliding  electron-po9itron  beams  cause  both  charge  and 
current  perturbations  giving  rise  to  mixed  polarizations.  Furthermore,  for 
the  case  of  colliding-beams  with  radial  dimension  smaller  than  the  colli¬ 
sionless  skin  depth  c/Wp,  the  finite  geometry  becomes  important  and  the 
usual  assumption  of  infinite,  homogeneous  medium  is  no  longer  valid.  In 
this  paper,  we  treat  the  f ilamentation  instability  of  colliding 
electron-positron  beams  with  finite  geometry  effects  included.  For  sim¬ 
plicity,  we  assume  in  this  section  that  this  colliding  beam  is  straight  and 
infinite  along  the  axial  direction. 

The  analysis  is  carried  out  within  the  framework  of  both  the  rigid 

beam  model  and  the  Vlasov -Maxwell  equations.  As  illustrated  in  Fig.  1,  the 

equilibrium  configuration  consists  of  intense  relativistic  electron  and 

positron  beams  propagating  opposite  to  each  other  with  axial  velocity 

Bee  for  the  positron  beam  and  Bee  for  the  electron  beam,  where 
p  ~z  r  e  ~z 

is  a  unit  vector  along  the  z-direction  and  c  is  the  speed  of  light  in 
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vacuo  and  0p— 6g.  Moreover,  both  beams  have  the  same  radius  R^  and  the 

o 

same  characteristic  energy  me  .  It  is  also  assumed  that  the  ratio  of 
the  beam  radius  to  the  collisionless  skin  depth  c/Wp  is  small,  i.e.. 


v  .  2  , 

- 1  =  N  —  ~  —  «  1 
Yb  j  me2  Yb 


where  j=e,p  denote  electrons  and  positrons,  respectively,  is  Budker's 

r°°  0 

parameter,  =  2it  dr  r  n.  (r)  is  the  number  of  particles  per  unit 
axial  length,  n^Cr)  *s  t*ie  equilibrium  particle  density  of  beam  component 
j,  -e  and  m  are  the  charge  and  rest-mass,  respectively,  of  electron.  As 
shown  in  Fig.  1,  we  Introduce  a  cylindrical  polar  coordinate  system 
(r,S,z).  All  equilibrium  properties  are  assumed  to  be  azimuthally  sym¬ 
metric  (3/39=0)  and  Independent  of  axial  coordinate  (3/3z=0). 

A.  Rigid  Beam  Model 

In  order  to  illustrate  the  physical  mechanism  of  this  f ilaraentation 
instability,  we  carry  out  the  stability  analysis  in  this  section  within  the 
framework  of  a  "rigid  beam"  model.  For  the  purpose  of  analytic  simplica- 
tion,  we  also  specialize  to  the  case  of  sharp-boundary  profiles  in  which 
the  equilibrium  density  profiles  are  rectangular,  i.e., 


"5<t)=  {o! 


n  =  const. 


0<r<R,  , 

D 


otherwise , 


where  j«e  and  p.  Making  use  of  Eq.  (2),  it  is  straightforward  to  show 
that  the  equilibrium  radial  electron  field  produced  by  particles  of  species 


j  is  given  by 
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{2ne  n  r  ,  0<r<R,  , 

3'J  2 

2fle  .n  ,R,/r  ,  r>R,  , 

J  J  b  b  * 


(3) 


where  ej  is  the  charge  of  particles  of  beam  component  j  (i.e.,  ej=-e  for 
j»e  and  ej»e  for  j-p) .  Similarly,  the  equilibrium  azimuthal  magnetic  field 
produced  by  particles  of  species  j  can  be  expressed  as 


r  271  VjV  ’ 

(r)  =  i 

^  2ire,n,b  .R^/r  , 
j  j  j  b 


J6 


0<r<Rb 


r>r 


b  » 


where  Vj*0jc  *s  the  axial  drift  velocity  and  c  is  the  speed 
vacuo . 

In  the  subsequent  analysis,  we  introduce  the  center  of 
coordinates  (Xj,Y^)  for  the  beam  component  of  specis  j.  In 
state,  we  assume  that 


(4) 


of  light  in 

mass 

the  equilibrium 


(Xj.Yj)  =  (0,  0)  ,  (5) 

for  j=e,  and  p.  It  is  also  assumed  that 

2  2  2 

Xj  +  Yj  «  R^  •  (6) 


The  restriction  to  small  perturbation  amplitudes  makes  the  subsequent  sta¬ 
bility  analysis  tractable.  The  transverse  motion  of  a  single  particle  of 
species  J  is  determined  approximately  from 


m , 


dt 


2  li 


c  dt  ~j 


BJ 


(7) 
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where  =  (x^.y^)  is  the  position  coordinate  for  a  particle  of  species  j 
and  E  and  g  are  the  total  electric  and  magnetic  fields,  and  m^Y^m  is  the 
relativistic  mass.  Assuming  E  and  B  can  be  approximated  by  their  equili¬ 
brium  values,  we  substitute  Eqs.  (3)  and  (4)  into  Eq .  (7).  The  equation  of 
motion  for  the  x  direction  can  be  expressed  as 

d2 

mj  ^2  xj  =  2"ej  £  Vk  1  -  VkAXj  "  V  ‘  (8) 

Neglecting  momentum  spread,  Eq .  (8)  can  be  averaged  over  the  beam  cross 
section.  After  some  straightforward  algebra,  we  obtain  the  approximate 
equation  for  average  motion  on  the  x  direction, 


Similarly,  the  equation  for  average  motion  in  the  y  direction  is  given 
by 


6  .b 

J  k 


(10) 


Defining 


Z.  =  X.  +  iYJ  (11) 

J  .)  j 


and  making  use  of  Eqs.  (9)  and  (10),  we  obtain 


dt 


2  Zj 


2ne 


BAn2J 


-  zk>  • 


(12) 


We  seek  oscillatory  wave  solutions  to  Eq .  (12)  of  the  form 


Zj  =  Zj  exp  iiikzi.z  +  i^c  tj  -  (Dtjj 


where  ui  is  the  complex  eigenf requency ,  Z^  =  const  is  the  perturbed  ampli¬ 
tude,  and  the  axial  wavenumber  kz  is  limited  to  the  range 


2  2 

k  R,  <  1  • 

z  b  ~ 


Equation  (14)  assures  the  approximate  validity  of  Eq .  (12)  for  wave  pertur 
bations  with  S/cJz^O.  Substituting  Eq.  (13)  into  Eq.  (12),  we  obtain 


„  2ire.  2ne. 

6  c)  M  VkU-6  BkJ,  Z  ■T‘|«A(WJV  \  - 

jk  J  jk  J 


Equation  (15)  gives  two  homogeneous  equations  relating  the  amplitudes  Z^ 

A  A 

and  Z  .  Setting  the  determinant  of  the  coefficients  of  Z,  equal  to  zero 
P  J 

gives  2X2  matrix  dispersion  equation  that  determines  the  complex  eigenfre- 
quency  w.  After  some  straightforward  algebra,  we  obtain  the  dispersion 
relation 


i(u>  +  k  cj^  -  a)  ^  !n  /n  )  j  if  a)  -  k  cj2-  j  =  0)^  [n  /n  j  , 
1  v  z  pe  p  e  v  z  pe  pe  pe 


1 2  2 


2  2" 

where  ft)  ■  4ire  ne/y^m  is  the  electron  plasma  frequency-squared  and  use 
has  been  made  of  8e=-6  =1,  which  is  consistent  with  present  experimental 
parameters . 

Assuming  that  both  electron  and  positron  beams  have  the  same  density. 


-  »•*  *»-■  T—..  *nw.  BfWH 


in*  l  Will 
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and  defining 


a  =  k  c  +  <jj  ,  h  =  kc-^)  , 

z  pe  z  pe 


we  simplify  the  dispersion  relation  in  Eq  .  (16)  as 


2  2  2  2  4 

i  <jJ  -a  i ,  ui  -b  ,  =  u>  , 

pe 


which  provides  a  necessary  and  sufficient  condition 


(17) 


(18) 


4  2,  2 

to  >  a  b 
pe 


,22 

=  i  k  c.  -uj 
z 


2  2 
pe' 


(19) 


for  Instability.  For  the  unstable  branch,  the  perturbation  is  purely 
growing  with  the  growth  rate 


w 

l 


7  9  2 

,  ta“-b  '  .  4 

Imo  =  |  l  ; — - — ,  +  w  j 


1/2 


2,2  1/2 

a  +b  i 


(20) 


The  maximum  growth  rate  of  instability  can  occur  at  a=0  or  b=0,  thereby 
giving 


1“,  J 

i  m 


=  15 


1/2  ,>1/2, 


pe  - 


0.  5  oo 


pe 


(21) 


For  colliding  beams  interacting  over  a  finite  distance  L,  the  axial 
wavenumber  kz  is  kz=2Tin/L  where  n=l,2,....  In  this  case,  the  condition  for 
a=*0  becomes  T»pe/c=2un.  The  finite  interaction  length  also  imposes  a 
severe  condition  for  the  instability  to  grow  significantly  before  the  beam 
exit.  Although  a  small  growth  of  perturbations  during  one  individual 


1 1 


interaction  of  electron  and  positron  beams,  we  expect  that  due  to  this 
f ilamentation  instability,  the  repating  interactions  between  both  beams 
eventually  convert  the  longitudinal  energy  of  the  beams  into  the  transverse 
energy  of  beams  and  the  field  energy  of  perturbations,  thereby  broadening 
the  beam  cross  section  and  leading  to  reduction  of  luminosity.  However, 
the  analysis  of  broadening  of  beam  cross  section  is  particularly  difficult 
and  is  currently  under  investigation  by  the  authors. 


B .  Vlasov  Description 


In  the  previous  section,  we  have  investigated  the  stability  properties 
of  dipole  oscillation  in  the  transverse  instability  for  the  electron  and 
positron  colliding  beams,  within  the  context  of  rigid  beam  model.  Although 
a  dipole  oscillation  in  a  rigid  beam  model  provides  a  simple  instructive 
description,  it  is  necessary  to  investigate  stability  properties  for 
perturbations  with  high  azimuthal  harmonic  number  i>2  within  the  framework 
of  the  Vlasov-Maxwel 1  equations. 

For  beams  of  well-defined  energy  and  momentum,  an  equilibrium 
associated  with  the  steady-state  (3/9t=0)  beam  distribution  function. 


f°(H,Pft,PJ  «  6[H-wP„— Y  mC2jd[  P  -Y  m3  ,  cj 


2^,  m 

o 


(22) 


is  particularly  suited  for  stability  analysis,  where  the  total  energy. 


u  ,'24  2  2,  ,  ,  , 

“  =  (.m  c  +  c£j  +  e^o(r)  , 


(23) 


the  canonical  angular  momentum, 
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pe  =  rPu  * 


(24) 


and  the  axial  canonical  momentum 


pz  +  ^j/cjAjr)  , 


(25) 


are  the  three  single-particle  constants  of  the  motion  in  the  equilibrium 

fields,  and  o>,  is  the  beam  rotational  frequency  of  species  j  and  Y.  is  a 
J  J 

constant.  In  Eqs.  (23)-(25),  $Q(r)  is  the  equilibrium  self-electric 

g 

potential,  A^(r)  is  the  axial  component  of  vector  potential  for  the 
azimuthal  self-magnetic  field,  and  g  =  (pr>PQ,Pz)  denotes  mechanical 
momentum  and  is  related  to  the  particle  velocity  v 
hy  l+E  /m  c  j 

Since  the  r-Q  kinetic  energy  of  particles  is  small  in  comparison  with 

O 

the  characteristic  energy  Ybmc  ,  it  is  straightforward  to  show  that  the 
term  H-WjPq  in  Eq .  (22)  can  be  approximated  by^ 


H  -  <y(.  =  yhmc  + 


+  ^7  \  a  r  » 


2Y,m  2  'b  j 


(26) 


where 


t;  -  C i -b p ) “ 1 1 


<1 


<>i  ■  Pf  +  lpo  ■  Ybl”jr) 


and 


'"jwj 


=  — U)  —  — 

j  Y,  in 


2ne  ^ 
b"‘  k 


(27) 


In  Eq.  (27),  the  laminar  rotation  frequency 


j 


is  defined  by 


u 

j 


2*e  . 

7^  J  nk  ek 

b  k 


ll_9AJ 


1  1/2 

J 


(28) 


Substituting  Eq .  (26)  into  Eq .  (22),  we  find  the  equilibrium  particle 
density  profile 


o ,  . 
(r) 


/d3p  f°(H,Pe,Pz) 


{ 


J 

0  , 


0<  r<R,  , 
b 

otherwise. 


(29) 


where  the  beam  radius  is  defined  by 


K  - 


(30) 


for  j=*e,p.  Equation(30)  ensures  that  the  electron  and  positron  beams  have 
the  common  beam  radius  R^.  It  is  important  to  note  from  Eqs.  (27)  and  (30) 
that  the  radially  confined  equilibrium  exists  only  for  the  rotational 
frequency  satisfying 


-U)  <  to .  <  U)  , 

j  J  j 


(31) 


Additional  equilibriium  properties  associated  with  the  distribution 
function  in  Eq.  (22)  are  discussed  in  Ref.  7. 

In  order  to  obtain  the  dispersion  relation  for  f ilamentation 


instability  of  the  electron  and  positron  beams,  we  make  use  of  the 
linearized  Vlasov-Maxwell  equations.  For  perturbations  with  azimuthal 


harmonic  number  l  and  axial  wavenumber  k^ ,  a  perturbed  quantity  <5 (x , t ) 

can  be  expressed  as  <S4>  (xjt)*!1  (r)exp{i(£b+k  z-wt)}  ,  where  w  is  the  complex 

eigenf requency •  The  present  stability  analysis  is  carried  out  in  long 

2  2  2 

parallel  wavelength  and  low  frequency  perturbation  satisfying  +1, 

2  2 

|wR^/c|  «l  +1  4  With  this  assmption,  the  axial  components  of  perturbed 
field  Ez(r)  and  Bz(r)  are  negligible  and  the  Maxwell  equations  of  perturbed 
potentials  can  be  expressed  as 


fi  3  a  i ,  4 

17 -~J  y(r)  = 

r 


-  4np(r) 


/i  a  3  4  4*  :  ,  , 

r  ■=■  -  —  j  A(r)  =  -  —  J  (r) 
vr  3r  3r  2'  c  z 


where  ^(r)  is  the  perturbed  electrostatic  potential,  P  (r)  is  the  perturbed 
charge  density,  A(r)  and  Jz(r)  ate  the  axial  components  of  the  perturbed 
vector  potential  and  current  density,  respectively.  Components  of 
perturbed  fields  can  be  expressed  in  terms  of  $(r)  and  A(r)  as 
E0  «  -U<Kr)/r,  Er(r)  =  -  (3/3r)i(r),  Br(r)  =  iU(r)/r,  and 
Bg(r)  =  -(3/3r)A(r). 

In  order  to  calculate  perturbed  charge  and  current  densities,  we  solve 
the  linearized  Vlasov  equation  to  obtain  the  perturbed  distribution 


function 


fjV.E) 


e.Y.m  3f 
J  b 


P,  3p 


3^  {‘'J 


(r)  +  (w-£w  -k  3  c j 
J  2  J 


dT  i  ^.(r  )exp[i  £  (0  -0)-i(o>-k  0  c)t]|  , 
J  2  J  J 


where  the  perturbed  electrostatic  potential  v^(r)  in  the  frame  of 

reference  moving  with  velocity  tfjC  is  defined  by  ^  (r)*^ (r)-d ^A(r)  and  use 

has  been  made  of  p  />,  m-e.  consistent  with  Eq .  (1).  It  is  useful  to 
z  h  j 

introduce  the  polar  momentum  variables  (p  ,?)  in  the  rotating  frame 

defined  by  p  +  r.mw.y  *  P,  cosy,  p  -  v.m^.x  =  p,  sin<y.  Note  also  that 
x  b  j  1  y  b  j  ri 

the  Cartesian  coordinates  (x,y)  are  related  to  the  polar  coordinates  (r,'J) 
by  x^r  cosd  and  y=r  sind.  In  this  context,  the  transverse  equation  of 
motion  of  particles  can  be  expressed  as^ 


x  (t  )  =  (  !/<*>  ;L(Pj^/^bml  cos<*>  sin  .>.r  -  ru^  sind  sinw^  t  +  cosb  cos^  jt  i 


y  (t)  “  (l/;»>.)f(p  /i,m)sinysin>u.T  +  ru  .cosjs  i  n  .  t  +  r.i .  s  i  nvcos .  . 

J  L  b  J  J  J  J  .1 


where  T»t  -t,  and  the  harmonic  frequency  w  is  defined  in  Eq.  (28). 

Upon  integration  of  Eq.  (34),  the  perturbed  charge  density  can  be 
found  to  be 


P(r)  -  2ne  \  y  m  dp,  p.  dp - L  [  V .  (r )  +  (w-tw  -k  d  cj  I  j  . 


1  1  1  Z  Pi  Pi  J 


j  *  i  '  y 


where  the  orbit  integral  Ij  is  defined  by 


I ,  -  i 


^  [  dT  4',  (r')exp  {  iL  £ iQ  -0)-(w-k  0  ,cjx]}  . 


Similarly  the  perturbed  axial  current  density  can  be  obtained.  For 
analytic  traceability,  we  will  consider  here  a  class  of  special  solutions 


J®  -mm*  « 
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for  which  the  perturbed  charge  and  current  density  are  localized  on  the 
beam  surface,  t.e.,  equal  to  zero  except  at  r«R^.  More  general 
perturbat lone .  particularly  the  body  wave  perturbations,  are  to  be 
presented  In  a  subsequent  publication.  In  this  case,  it  follows  from  Eqs. 
(32)  and  (33)  that  the  function  ^j(r)  has  the  simple  form  Vj  (r)  = 

$(r)-b  A(r)  *  C^r*  for  0<r<R},*  Substituting  Eq .  (35)  into  Eq .  (37)  it  is 

readilv  shown  that 

i  1  .  (  r )  [0 

l.  =  --L- -  a  i  exp  |  -  i  (.-k  :•  .<);  J  [(..+<„.  )exp(  iw .  t  ) 

i  ,  >  1  z  i  .1 

U‘V  '  (38) 

-  (  .  .  -  .  .  )  exp  (  -  !  ■•  .  '  )  i 
1  I  i 

After  some  straightforward  algebra  that  utilizes  Eqs.  (22),  (36)  and  (38), 
Eq.  (.321  can  he  expressed  us 


1  f  r  f  -  fI  V  (D  - 
-  3r  ?r  t 


l  *  (r)  --Pj  r  (u)  6(r-Rb)  , 

i  i  ^ .  K  J 

i  b 


2  2  ~ 

where  u)  »  4ne  n./y,m  is  the  plasma  frequency  -  squared  of  beam 
PJ  1  b 

9 

component  j,  »  (u<  ^  )  ('o  ^  -v  ^  )  is  defined  in  Eq  .  (27)  and  r  j  (co )  is 

defined  by 


r  U) 


:)  ^  i,  U'.-k  r-.C  r  U!  .  +V  .  )  ° 

=  -  1  +  |  J  .  J  V  ‘  !  .  .  .  _J  7  J  .....  J  J 

(.  2uj  j  n  !  (  f  -n)  !  (*>-k  !•  .  c+ » <i  .  -  2n.-.( .  Iuj.-ll' 


[  2ui .  j 
v  j  '  n=0 


z  1  .) 


Similarly,  Eq.  (33)  can  be  expressed  as 


13  3  e‘ 


r  -3-  -  - -j]  A(r)  =  -  l  P1«.(r)  f .  (u>)6  (r-R^)  . 


where  use  has  been  again  made  of  the  approximation  p,/y  m  —  &jC  consistent 


with  Eq.  (1) 
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As  the  right-hand  sides  of  the  coupled  differential  equations  (39)  and 
(41)  are  equal  to  zero  except  at  the  surface  of  the  beam  r«Rjj»  they  can  be 
solved  In  a  straightforward  manner  to  give 


C 


k 


)'d 


W 


r.(u,)c.  . 

i  ) 


(42) 


In  the  case  when  the  beams  are  located  Inside  the  cylindrical  conducting 

wall  with  radius  Rc,  the  term  In  the  left-hand  side  of  Eq.  (42)  Is 

replaced  by  [l-^jj/Rp)^]-^^.  Note  that  the  absolute  value  of 
2  2 

(o>)/il  in  Eq.  (42)  is  of  the  order  of  unity  or  less.  It  follows  from 
Eq.  (42)  that  the  condition  for  a  nontrivial  solution  (Cj  not  all  zero)  is 
given  by 


1  -  <“pp  V/t2  “p  rp<“)  re‘“>  '  °  • 


(43) 


—2 

where  use  has  been  made  of  3  --3  =1  and  Y,  <<1,  which  Is  consistent  with 

p  e  b 

present  experimental  parameters.  Equation  (42),  when  combined  with  Eq. 
(40),  constitutes  one  of  the  main  results  of  this  paper  and  can  be  used  to 
Investigate  f ilamentation  stability  properties  for  a  broad  range  of  system 
parameters. 

As  an  example,  we  restrict  the  investigation  of  dispersion  relation 
(43)  to  the  case,  where  both  beams  are  In  a  cold  fluid  rotational 
equilibrium  characterized  by  *  A  careful  examination  of  expression 

for  r j  (oj)  show  that^ 


lim  , 

CO  *>±U) 

j  j 


vl 


2in 


2  V“> 


121 


(44) 


2(o)-k  B.c+tu). )  [at-k  8 ,c+(£-2)w.  ] 
z  J  J  z  J  J 


Therefore,  in  a  cold  fluid  limit,  the  dispersion  relation  in  Eq.  (43)  can 
be  considerably  simplified-  After  some  algebraic  manipulation  we  can  show 
that  fof  the  fundamental  mode  perturbation  (i.e-,  *.=1),  the  dispersion 
relation  in  Eq.  (43)  is  identical  to  Eq .  (16)  obtained  within  the  framework 
of  rigid  beam  model.  The  stability  analysis  of  Eq .  (43)  for  a  broad  range 
of  harmonic  number  t  and  rotational  frequency  is  currently  under 
investigation  by  the  authors.  Nonlinearly  the  beams  become  filamentated 
first,  then  the  current  filaments  of  the  same  sign  attract  each  other  to 
form  a  broader  beam.  Finally,  we  conclude  this  section  by  pointing  out 
that  the  understanding  in  broadening  in  beam  cross  section  by  repeating 
interactions  of  beams  is  not  completed  yet.  And  this  area  is  currently 
under  investigation  by  the  authors. 
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I [ I •  COLLECTIVE  SELF-FIELD  EFFECTS 

In  this  section,  we  examine  the  influence  of  the  collective  self- 

fields^  on  the  electron  and  positron  colliding  beams  In  the  storage  ring. 

While  the  forces  of  the  self-generated  electric  and  magnetic  field  of  a 

highly  relativistic  electron  (positron)  beam  on  an  electron  (positron) 

-2  ,  -2 

cancel  out  to  order  0(y  ),  i.e.,  F,  +  d  BAaO('f  )  the  forces  of  the 

r  c  u 

electric  and  magnetic  fields  of  the  electron  beam  on  the  colliding 
positrons  are  additive  leading  to  radial  acceleration.  This  effect  of  the 
collective  self-fields  of  one  species  on  the  other  species  of  the  colliding 
beams  Imparts  considerable  transverse  energy,  thereby  substantially 
increasing  the  beam  transverse  dimensions  upon  collision.  In  order  to  make 
the  problem  simple,  we  assume  that  the  colliding  section  of  the  storage 
ring  is  straight.  The  theoretical  analysis  is  carried  out,  distinguishing 
the  two  cases,  where  (a)  the  particle  motions  are  in  a  very  coherent  orbit 
and  (b)  the  randomness  dominates  the  operational  condition  of  storage  ring 
(e.g.,  Incoherent  collision  location  by  fluctuation,  etc.).  In  either 
case,  it  is  found  that  the  self-fields  effects  play  a  dominant  role  in  the 
stability  behavior  of  transverse  orbit  or  the  expansion  of  beam  cross 
section.  For  present  experimental  parameters^  at  DESY,  the  cross  section 
of  the  beam  can  be  expanded  to  ten  times  of  its  original  area  within  5 
milliseconds  operational  time.  Without  loss  of  generality,  we  assume  in 
Fig.  1  that  the  front  edges  of  both  beams  arrive  in  z=0  at  time  t*>0. 

The  axial  orbit  of  particles  of  beam  component  j  is  given  by 

(45) 

I  1 

where  the  initial  position  /,  is  rest  rioted  to  satisfy 


*,(/  +■  ,1.1  •  0  . 

1  I  1 


(4  b) 
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Here  £j  *  sgn  t*.  and  e.  is  the  charge  of  the  particles  of  beam  component  j. 
The  particle  density  profile  of  beam  component  j  is  expressed  as 

n°(r,z.t)  =  n. (r,z)u[(&.rt-z) (z+e.l-BjCt) )  ,  (A  7) 

where  the  Heaviside  step  function  U(x)  is  defined  by 

(0.  x  <  0  , 

U(x)  =  <  (48) 

(  1,  x  >  0  . 

For  a  specific  choice  of  the  beam  density  n.(r.z)  in  Eq .  (47),  the 
potentials  for  the  self-fields  are  to  he  calculated  from  the  Maxwell 
equations.  The  Poisson  equation  can  be  approximated  bv 

18  8 

r  Jr  r  Jr  4l(r'z’t)  =  ~  4n^e  n^r.z.f)  ,  (49) 

•  .1  J 


where  ()>(r,z,t)  is  the  self-electric  potential.  In  obtaining  Eq.  (49),  we 

2  2 

neglect  the  term  proportional  to  3  <Ji/3  z,  under  the  assumption  that  the 
axial  length  L  of  the  beam  is  much  larger  than  the  beam  radius  and  the 
effects  of  the  leading  edge  of  the  beams  are  thus  neglected.  Furthermore, 

the  z-component  of  the  VxBs(x)  Maxwell  equation  is  expressed  as 

a.  -v  ^ 


1  J_  J 

r  3r  r  3r 


A^r.z.t) 


4irye.t5.ri?(r,z,t)  , 

J  1  1  1 


(50) 


where  Az(r,z,t)  is  the  z-component  of  the  self-vector  potential.  Other  com¬ 
ponents  of  the  vector  potential  are  negligible  because  of  Eq.  (1).  Defining 

the  effective  self-potential  (r , z, z . )  =  <f>-6 .As,  and  making  use  of  Eqs.  (45) 

J  1  1  2 

(47),  (49),  and  (50),  we  have 


_3_ 

3r 


<|/!(r,z,z  )  =  -  8 ire  -- 
J  1  k  r 


r 

dr'r'nk(r' ,z) 


x  U[(zj-2z)(2z-z.+t.kL)l  , 


(51) 


i 
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where  k# j .  In  obtaining  Kq .  (51),  use  lias  been  made  of  y  ~2  =  (i_g2)<<} 

b  p 

In  order  to  make  the  problem  simple,  we  carry  out  the  analysis  in  the 
average  applied  field  provided  externally  by  the  periodic  quadrupole  mag¬ 
netic  field,  similar  to  that  used  in  the  previous  study*\  In  this  regard, 
the  applied  focussing  force  can  be  obtained  from  the  axial  component  of 
the  effective  vector  potential 

Az  (r)  =  -(Tfbm/2eBp)  wf  r  (52) 


where  is  the  focussing  oscillation  frequency  determined  by  the 
quadrupole  field  gradient. 

The  total  energy  of  particles  oi  the  beam  component  j  is  given  by 
24  22 

H  =  (me  +c  p'")  +  e  ,6(r,z,t),  (53) 

where  the  lower  case  p  denotes  mechanical  momentum  and  is  related  to  the 

2  22  1/2 

particle  velocity  v  by  v  = p/m(l+p  /m  e  )  .  Since  the  r-9  kinetic 

'V  X,  'X,  'V 

energy  of  particles  is  small  in  comparison  with  the  characteristic  energy 
2 

■y^mc  and  Vj/y^<<l  in  Eq.  (1),  it  is  straightforward  to  show  that  Eq .  (53) 
can  be  approximated  by 


2 

2  px+py  s  1^2 

H  =  Vc  +  27V  +  ei+j(r’Z’Zj)  +  2  Ybm  wf  r“  ’ 
b 


(54) 


2  2 

where  y^  =  C 1— 6 ^ )  .  From  Eq.  (54),  we  obtain  the  equation  of  motion  for 


Z(t)  =  x(t)  +  iy(t) 


(55) 


1/2  2 
where  i=  (-1)  •  Making  use  of  Eqs .  (45)  and  (55),  and  6^-1,  the 

equation  of  motion  for  particles  of  the  beam  component  j  is  given  by 


g  2  [Y  U)f 

— |  +  -H - 2  *  ~2  dr,r,nk(r’ ,z)U[  (7.  ,-2z)  (2z-z  +ekL)  ]  +  —  z  *  0  (56) 

dzZ  y,mc  r  'o  1  ’  c 

D 
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where  k^j  and  z.  is  defined  in  Eq.  (45).  Equation  (56)  determines  the  transverse 
position  of  particles  of  beam  component  j,  thereby  providing  the  information 
of  the  particle  density  n .  ( r ,  /, )  ,  which  in  turn  governs  the  equations  of 
motion  for  particles  of  the  beam  component  k.  In  this  regard,  the  coupled 
differential  equation  (56)  for  j=e  and  p  can  be  used  to  investigate  the 
temporal  profile  evolutions  of  various  beam  properties  for  a  broad  range  of 
initial  parameters. 

As  an  example,  we  consider  a  tenuous  positron  beam  satisfying 

(O2  <*  (c/l.r  ,  (57) 

pp 

-2  -  2 

where  =  4nn^e  is  the  average  positron  plasma  frequency-squared 

in  the  laboratory  frame.  Equation  (57)  assures  that  all  the  electrons 
move  on  the  straight  paths  with  constant  radius  r  during  the  collision. 

Assuming  the  electron  density  profile  as 


,  ,  fv  r<Rb  • 
".('.a  ■  1 

,  otherwise. 


the  transverse  equation  of  motion  for  positron  can  be  expressed  as 

,  2  2 
d  2?  W 

— j  +  -jp  Z  U[(z  -2z)(2z-z  -L) ]  +  Z  =  0  , 
dz  c^  p  p  c“ 


2  "2 

where  u>  =  4nnee  /-y^m  is  the  electron  plasma  frequency-squared. 

Without  loss  of  generality,  we  assume  that  there  is  one  pair  of  electron 
and  positron  beams  in  the  entire  system,  thereby  indicating  that  the  whole 
storage  ring  can  be  represented  by  two  focusing  sectors.  Each  sector  consists 
of  a  self-beam  focusing  set  (the  region  in  which  beams  collide)  and  an 
applied  focusing  set.  The  subsequent  analysis  is  carried  out  distinguishing 
the  two  cases:  (a)  the  positrons  move  on  a  very  coherent  orbit,  and  (b)  the 


axial  location  of  collision  as  well  as  the  beam  length  fluctuates  incoherently, 
thereby  the  ensemble  average  can  be  feasible. 


23 


A -  Stability  Analysis  of  Coherent  Positron  Orbit 


The  stability  properties  of  individual  particle  orbit  can  be  determined 
from  the  transformation  matrix  of  one  sector  for  a  very  coherent  positron 
orbit.  Assuming  that  a  positron  has  an  initial  condition  Z=Z^  and 
Z'  =  (dZ/dz)  =  Zj  at  z = z^/2,  it  can  be  shown  from  Eq .  (59)  that  the  transverse 
orbit  of  this  positron  is  given  by 

Z  =  Z  cos  [  (wT/c)  (z-z  /2)  ]  +  (Z  c/m  )sin[  (0jT/c)  (z-z  /2)  ]  ,  (60) 

IT  p  II  1  p 

2  2  1/2 

for  z  /2<z<z  /2+L.  Here  the  frequency  =  (w  +wr)  .  From  Eq.  (60) 
p  n  1  T  pe  f 

it  is  also  straightforward  to  show  that  the  transverse  position  Z^  and  orbit 
slope  of  positron,  when  it  emerges  from  the  right-hand  side  of  the 
electron  beam,  is  given  by 

(Z2\  /  cos(aiTL/2c)  (c/wT)sin(wTL/2c)\  /  zA 

,Z2  /  \-(uT/c)sin(uT  /2c)  cos  (m.^L/Zc )  /  \  Z^  / 


Similarly,  when  the  applied  focusing  section  has  been  traversed,  the 
position  and  orbit  slope  are  given  by 


/  Z  3  \  /  cos<f> 

\Z  \-(uf/c)i 


(c/mf) sint\  /Z 


cos# /  \Z 


where  the  phase  shift  <f>  =  u>f (S-L) /2c  and  S  is  the  length  of  the  whole  circum¬ 
ference  of  storage  ring. 

Therefore,  from  Eqs.  (61)  and  (62),  we  obtain  the  trace  of  the  trans— 


I 


formation  matrix  M  for  a  sector 
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Tr  M  =2  cos'!’  cost  ; — )  -  (  ~  +  ~  )sin‘I>  sin(----j  ,  (63) 

ZC  Ul  - 

which  is  the  sum  of  the  elements  of  the  principal  diagonal  of  the  transforma¬ 
tion  matrix  M.  The  necessary  and  sufficient  condition  for  stable  transverse 

.  .  •  8 
orbit  is 


,1  ,, 


r  M  .  '  I 


(64) 


As  a  typical  example  in  the  present  experiment,  we  consider  the  system 
9  7 

parameters  w  £  =  10  rad/sec,  i,  =  2cm,  and  _r-  2*10  rad/sec.  Substituting 

these  parameters  into  Eq.  (63)  gives  approximately  Tr  M/2  -  cos!  -sin!, 

which  violates  the  inequal  it  v  in  F<(.  (64)  for  the  range  (n-0. 5)tt  <  f  <  n” , 

where  n  is  an  integer.  We  therefore  conclude  that  the  collective  self-fields 

effects  (to  )  of  the  electron  and  positron  colliding  beams  play  a  signifi- 
pe 

cant  role  in  the  stability  behavior  of  transverse  particle  orbit. 


B.  Expansion  of  Beam  Cross  Section  witli  Ensemble  Average 

In  order  to  investigate  die  expansion  of  beam  cross  section  for  uncon¬ 
trollable  collision  (incoherent  collision  location,  ere.),  we  define 

2  *  ?  »  ?  sfc 

rI  =  Z1Z1  +  <c/“f )"Z1Z1  >  (65) 

which  represents  the  maximum  radial  deviation  from  the  axis  of  symmetry 
before  collision.  In  Eq. (63)  ,  the  asterisk  (*)  denotes  the  complex 
conjugate.  During  the  collision  (z  /2  <  ?■  <  z  /2  +  I.) ,  the  transverse  orbit 
of  a  positron  in  Eq.  (60)  can  also  be  expressed  as 

Z  =  A  cosf  (u>T/c)  (z-z  /2)  +  «]  ,  (66) 

where  A  is  the  maximum  amplitude  and  i  is  the  initial  phase  angle  which  is 
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defined  by  a  =  tan  ^  [ - (c/m^, )Z ^/Z ^  j .  The  maximum  radial  deviation  for  ranee  z 
satisfying  z  /2<z<z  /2+1,  is  determir.  _  from 

AA*  =  r  /[l+(iu  /to.)  ^sin^a )  ,  (67) 

J  pe  t 

2  1/2 

where  use  has  been  made  of  Eqs.  (65)  and  (66),  and  =  (u  +w^)  .  From 

I 

Eq.  (66),  the  positron  position  7,^  and  orbit  slope  can  be  expressed  as 


7-2  ~  A  cos  [  (ui,|,L/2c)  +  a  ]  , 

Z2  =  -A(wT/c)sin[  (w,_,L/2c)  +  a]  , 


thereby  giving  the  relation 


r2\2  l+(w  /wf)  sin  [ (mTL/2c)+a ] 


(W - 


2  2 

l+(w  /to,)  sin  a 
pe  t 


(68) 


from  which  the  maximum  radial  deviation  r^  after  collision  is  determined. 

Depending  on  the  phase  angle  a,  positrons  gain  (or  lose)  the 
transverse  energy  by  the  collision  according  to  r2/r^>l  (or  r2/r^<l).  The 
net  gain  of  the  transverse  energy  (or  temperature)  by  the  collision  is 
determined  from  the  phase  angle  average  of  Eq.  (68).  We  therefore  define 


9  9  .  r  2rr  l+(m  Ao  c)  ^sin^  [  (m  L/2c)-fn] 

-  2 ,  2  1  pe  f  T 

<r2/ri>=  2tt  da 

o 


2  2 

1+0.)  /uO  sin  a 
pe  i 


f  (a) 


(69) 


for  future  notational  convenience.  In  Eq.  (69),  the  phase  angle  distribution 

f(a)  is  a  positive  definite  function  normalized  bv  f  da  f (a)  =  2tt  For 

Jo 

uniform  distribution  (f=l),  we  obtain 


2  2 
1 4to  /2io  - 


2  0  (Urri1-*  /z.0). 

<r2/r>-  : » +;— f1-"8* 


wtL 


(1-ho  /to  ) 
pe  f 


)]• 


(70) 


2  2 

Evidently,  we  note  from  Eq.  (70)  that  the  value  tr^/r^  approaches  unity 


I 
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when  the  beam  length  (L)  or  density  (w  )  decreases  to  zero.  Moreover,  the 
2  2 

value  is  always  greater  than  unity. 

As  a  typical  example  in  the  present  experiment,  we  evaluate  Eq. (70) 

9  7 

for  w  =10  rad/sec,  1.  =  2cm  and  w-  =  2><10  rad/sec.  Substituting  these 
pe  t 

2  2 

parameters  into  Eq.  (70),  we  find  <r9/r^>  =  1. 025.  Therefore,  in  these 
particular  parameters,  the  cross  section  of  the  beam  is  increased  by  2.5 
percent  of  its  original  area  after  each  collision.  However,  we  assume 
that  the  positrons  are  uniformly  distributed  in  the  phase  angle  a  whenever 
beams  start  collision,  which  is  consistent  with  the  ensemble  average 

scheme.  The  cross  section  of  the  positron  beam  can  be  expanded  to  ten 

2  2 

times  of  its  original  area  for  c^/r^  =  1-025  after  100  times  collisions, 
which  corresponds  to  the  operational  time  (S/2c) logl. 025 =  5  milliseconds 
for  the  circumferential  legnth  S =  3xl0^cm  of  storage  ring. 
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IV.  CONCLUSIONS 

In  this  paper,  we  have  examined  the  f ilamentation  instability  and  the 
influence  of  the  collective  self-fields  on  the  electron-positron  colliding 
beams  in  the  storage  ring.  In  Sec.  II,  we  have  investigated  the  stability 
properties  of  f ilamentat ion  instability  of  electron-positron  colliding 
beam.  An  important  conclusion  of  this  stability  analysis  is  that  the 
typical  growth  rate  of  the  f ilamentation  instability  is  order  of  the 
electron  plasma  frequency,  thereby  severely  limiting  the  electron  density 
in  a  storage  ring.  Influence  of  collective  self-field  effects  on  the 
electron  and  positron  colliding  beams  has  been  investigated  in  Sec.  III. 
The  theoretical  analysis  has  been  carried  out,  distinguishing  the  two 
cases,  where  (a)  the  particle  motions  are  in  a  very  coherent  orbit  and  (b) 
the  randomness  dominates  the  operational  condition  of  storage  ring  (e.g.. 


incoherent  collision  location  by  fluctuation,  etc.).  In  either  case,  it 
has  been  found  that  the  self-fields  effects  play  a  dominant  role  In  the 
stability  behavior  of  transverse  orbit  and  the  expansion  of  beam  cross 
section. 
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FIGURE  CAPTION 


Fig.  1  System  configuration  and  coordinate  system. 
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